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1.  INTRODUCTION 

The  effort  to  be  described  in  this  report  was  performed  as  partial  fulfillment  of  two 
primary  objectives:  (1)  perform  computer  simulations  of  charged  particle  transport, 
energy  and  charge  deposition  in  satellite-borne  instrumentation  used  in  research  efforts  of 
the  Air  Force  Research  Laboratory/Space  Weather  Center  of  Excellence  (AFRL/VSBXR) 
to  detect  and  characterize  (by  type,  energy,  intensity,  etc.)  particles  associated  with 
ionizing  radiation  in  space;  (2)  transfer  this  simulation  capability  to  AFRL/VSBXR  and 
provide  advice  to  Air  Force  researchers  on  its  use;  and  (3)  perform  studies  of  transport 
code  predictions  of  grazing  angle  proton  scattering.  In  the  following  sections,  we  provide 
descriptions  and  examples  of  particle  transport  simulations  and  their  application  to 
problems  of  interest  to  AFRL/VSBXR. 

2.  ELECTRON  TRANSPORT  MODELING 

2.1  CRRES  Dosimeter  Models  -  MCNPX  and  ITS/ACCEPT 

A  series  of  flux  and  dose  response  functions  for  electrons  and  protons  were 
calculated  by  Auehampaugh  and  Cayton  [1]  for  electron,  secondary  photon  and  proton 
transport  in  the  CRRES  (Combined  Release  and  Radiation  Effects  Satellite)  dosimeters. 
The  calculations  were  made  in  1992  using  the  then  current  versions  of  MCNP[2]  for 
electron  and  photon  transport  and  LAHET[3]  for  the  proton  transport.  Since  that  time, 
the  MCNP  and  LAHET  codes  have  been  superceded  by  MCNPX[4].  Also,  our 
experience  has  shown  that  the  electron  transport  calculations  can  in  general  be  performed 
more  efficiently  using  ITS/ACCEPT[5]  rather  than  MCNPX.  For  this  reason,  as  well  as 
recognition  of  significantly  (since  1992)  reduced  computation  cost,  we  undertook  to 
repeat  some  of  the  electron/photon  transport  calculations  of  [1]  using  both  the  MCNPX 
and  ACCEPT  codes,  thus  providing  an  “in-house”  Air  Force  capability  to  perform  these 
and  other  simulations  as  may  be  needed  in  the  future. 

The  original  MCNP  geometry  files  from  [1]  for  the  CRRES_D1,  D2,  D3  and  D4 
dosimeters  were  supplied  by  D.  Brautigam  (AFRL)[6].  The  D1  and  D2  files  are  shown 
in  Figures  1  and  2,  respectively. 

We  used  the  visual  editor  program  VISED[7]for  MCNP  input  files  to  assist  in  the 
interpretation  of  the  geometry  files  shown  in  Figures  1  and  2.  This  visualization  is  shown 
in  Figures  3  and  4,  which  are  the  diagrams  for  the  CRRES_D1  cell  and  surface 
configurations,  respectively.  VISED  provides  the  color  rendering  to  distinguish  among 
the  five  materials  specified  by  the  “mX”  (X=l-Si;  2-A1;  3-Al203;  4-Ni;  5-W,Fe  alloy), 
records  in  the  geometry  files. 
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Geometry  for  CRRES  dosimeter  dome  1 

11  pz  0.00000 

-11  Feb  92 

12  pz -0.07620 

1 

1  -2.3311  -10-3 

13  pz -0.10668 

2 

1  -2.3311  -10  3  23-24  25-26 

14  pz -0.15748 

3 

011  -10  (-23:24:-25:26)  -2 

15  pz -0.31496 

4 

010-2 

16  pz -0.94996 

5 

2-2.700  2-1  10 

17  pz -1.58496 

6 

2-2.700  2-1  11  -10 

18  so  1.46812 

7 

012-11  -22 

19  S0 1.32000 

8 

013-12-5 

20  SO  50.0 

9 

3-3.70012-11  22-8 

21  pz  -0.0000001 

10 

2-2.70012-11  8-7 

22  cz  0.06350 

11 

2-2.70013-12  5-7 

23  px  -0.09525 

12 

4-8.90014-13-5 

24  px  0.09525 

13 

2-2.70014-13  5-7 

25  py  -0.09525 

14 

2-2.70015-14-7 

26  py  0.09525 

15 

5-18.30  17-15-6 

1  16 

2-2.700  17-156-7 

mode  p  e 

17 

2-2.700  17-167-9 

imp:p,e  1  1 9r  0 1 

18 

0  16-11  7-9 

ml  140001. 

19 

0  11  1  -18 

m2  130001. 

20 

0  ((11  18):(-11  9):-17)  -20 

m3  8000  -0.47075  1 3000  -0.52925 

21 

0  20 

m4  280001. 

22 

0  21 

m5  74000  -0.95  26000  -0.05 
sdef  sur=19  ccc=22  nrm=-1.  dir=d1 

1 

SO  1.22936 

par=3  erg=5.0 

2 

SO  1.01981 

sbl  -21  1 

3 

cz  0.05093 

spl  -21  1 

4 

cz  0.76200 

f108:p,e  1 

5 

cz  0.32385 

el  08  0.  .049  .051  .125  .193  .263  .336 

6 

cz  1.1 1125 

.408  .480  .549  .622  .694  .765 

7 

cz  1.22936 

.839  .910  0.979  1.020  1.28  1.91  2.56 

8 

cz  0.50800 

3.20  3.84  4.48  5.12  5.74 

9 

CZ  1.46812 

6.37  7.01  7.67  8.30  8.94  9.57  10.1 

10 

pz  0.04030 

41 .2  104.  t 
phys:e  25.5 

Figure  1.  CRRES_D1  Dosimeter  run  file,  LANL  Version[l],  for  MCNP[2]. 
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Geometry  for  CRRES  dosimeter  dome  2 
-11  Feb  92 

1  1-2.3311-10-3 

2  1  -2.3311  -10  3  23-24  25-26 

3  011  -10  (-23:24: -25:26)  -2 

4  010-2 

5  2-2.700  2-1  10 

6  2-2.7002-1  11-10 

7  012-11-22 

8  013-12-5 

9  3-3.70012-1122-8 

10  2-2.700  12-11  8-7 

11  2-2.70013-12  5-7 

12  4-8.900  14-13-5 

13  2-2.70014-13  5-7 

14  2-2.700  15-14-7 

15  5-18.3017-15-6 

16  2-2.700  17-15  6-7 

17  2-2.700  17-16  7-9 

18  016-117-9 

19  0111-18 

20  0  ((11  1 8):(-1 1  9):-17)-20 

21  0  20 
22  0  21 

1  so  1.74498 

2  so  1.15443 

3  cz  0.12741 

4  cz  0.76200 

5  cz  0.32385 

6  cz  1.54813 

7  cz  1.74498 

8  cz  1.14300 

9  cz  1.90500 

10  pz  0.04340 

11  pz  0.00000 

12  pz -0.07620 


13  pz -0.10668 

14  pz -0.15748 

15  pz -0.31 496 

16  pz -0.94996 

17  pz -1.58496 

18  so  1.90500 

19  S0 1.76000 

20  so  50.0 

21  pz -0.0000001 

22  cz  0.12700 

23  px -0.19685 

24  px  0.19685 

25  py -0.19685 

26  py  0.19685 

mode  p  e 
imp:p,e  1  1 9r  0  1 
ml  140001. 
m2  130001. 

m3  8000  -0.47075  1 3000  -0.52925 
m4  280001. 

m5  74000  -0.95  26000  -0.05 
sdef  sur=19  ccc=22  nrm=-1.  dir=d1 
par=3  erg=2.00 
sbl  -21  1 
spl  -21  1 
f108:p,e  1 

el  08  0.  .051  .066  .134  .217  .287  .356 

.434  .507  .579  .655  .735  .807 

.885  .955  1.035  1.26  1.95  2.58  3.23 
3.87  4.51  5.13  5.75 

6.38  7.02  7.66  8.28  8.91  9.83  10.21 
42.1  104.  t 
phys:e  5.5 


Figure  2.  CRRES_D2  Dosimeter  run  file,  LANL  Version  [l],for  MCNP[2]. 
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Figure  3.  VISED  [7]  rendering  of  MCNP  cell  geometry  configuration  [1] 
for  CRRES_D1,  corresponding  to  the  listing  shown  in  Figure  1. 


Figure  4.  VISED  [7]  rendering  of  MCNP  surface  geometry  configuration 
[1]  for  CRRES_D1,  corresponding  to  the  listing  shown  in  Figure  1. 
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2.1.1  CRRES_D1  MCNPX  Calculations 

The  authors[l]  used  a  cosine- weighted  source  (as  specified  by  the  “sdef,  sbl  and 
spl”  records  of  Figures  1  and  2).  A  cosine- weighted  current  source  corresponds  to  an 
isotropic  flux  of  electrons  incident  on  die  dome  surface.  In  general,  the  relationship 
between  current  j  and  flux  (j)  is 

j(jU)  =  JMp{n)  (1) 

where  ju  is  the  cosine  of  the  polar  angle  between  the  particle  direction  and  the  normal  to 
the  surface  of  incidence. 

In  the  LANL  calculations!!]  it  was  assumed  that  the  source  particles  originated  on 
a  hemispherical  surface  concentric  with,  but  slightly  larger  than,  the  dosimeter  dome. 
The  CRRES_D1  input  file  (Figure  1)  was  used  to  perform  energy  deposition  calculations 
in  the  silicon  dosimeter  (cell  #1  as  defined  in  Figure  1)  with  MCNPX  for  four  electron 
source  energies,  2,  5,  10,  25  MeV,  and  die  results  were  compared  with  those  obtained  in 
[1].  These  comparisons  are  shown  in  Table  1.  We  computed  the  energy  deposition  using 
two  methods:  1)  the  energy  deposition  tally  supplied  by  MCNPX;  and  2)  the  method  used 
in  [1],  a  summation  over  the  pulse  height  distribution  of  energy  deposited.  The  energy 
deposition  tally  capability  had  not  yet  been  incorporated  into  MCNP  when  the 
calculations  of  [1]  were  made.  Only  the  pulse  height  distribution  summation  method  was 
available  at  that  time. 
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Table  1.  Energy  Deposition  (MeV)  in  CRRES_D1  Dosimeter  (Cell  #1,  Figures  1,3) 


SOURCE 

ENERGY 

(MEV) 

MCNPX 
METHOD  1 
COSINE 
WEIGHTED 

MCNPX 
METHOD  2 
COSINE 
WEIGHTED 

ACCEPT 
METHOD  1 
COSINE 
WEIGHTED 

ACCEPT 
METHOD  2 
COSINE 
WEIGHTED 

AUCHAM- 
PAUGH  & 
CAYTON 

m' 

ACCEPT 
METHOD  1 
ISOTROPIC 

IN  COSINE 

ACCEPT 
METHOD  2 
ISOTROPIC 

IN  COSINE 

0.2 

2.232e-8 

8.79  le-9 

5.90e-10 

0 

0 

0.3 

1.108e-8 

0 

1.4  le-9 

9.512e-9 

0 

0.5 

2.698e-8 

8.79  le-9 

3.8  le-9 

1.708e-8 

8.79  le-9 

1 

4.445e-8 

2.258e-8 

2.33e-7 

6.560e-8 

2.468e-8 

1.5 

6.059e-5 

6.017e-5 

4.27e-5 

4.615e-5 

4.58  le-5 

2 

1.790e-4 

1.856e-4 

1.975e-4 

1.984e-4 

1.07e-3 

1.627e-4 

1.610e-4 

3 

3.114e-4 

3.106e-4 

1.86e-4 

2.580e-4 

2.558e-4 

4 

3.060e-4 

3.054e-4 

2.09e-4 

2.477e-4 

2.459e-4 

5 

2.898e-4 

3.193e-4 

2.668e-4 

2.73 2e-4 

2.56e-4 

2.166e-4 

2.185e-4 

6 

2.604e-4 

2.678e-4 

2.73e-4 

1.95  le-4 

2.009e-4 

r _ 

2.454e-4 

2.564e-4 

i 

2.73e-4 

1.805e-4 

1.86  le-4 

■ 

2.467e-4 

2.552e-4 

2.50e-4 

1.829e-4 

1.893e-4 

10 

2.828e-4 

2.846e-4 

2.380e-4 

2.899e-4 

2.69e-4 

1.679e-4 

1.737e-4 

. 

12,5 

2.324e-4 

2.41 5e-4 

2.77e-4 

1.570e-4 

1.689e-4 

15 

2.26  le-4 

2.43  le-4 

2.73e-4 

1.572e-4 

1.676e-4 

20 

2.156e-4 

2.372e-4 

2.63e-4 

1.523e-4 

1.583e-4 

25 

2.729e-4 

3.222e-4 

2.215e-4 

2.433e-4 

2.60e-4 

1.497e-4 

1.648e-4 

l 


Ref.  1  (Energy  Deposition)  /  Q 
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2.1.2  CERES_D1  ITS/ACCEPT  Calculations 

The  CRRES_D1  geometry  and  materials  data  given  in  Figure  1,  along  with  the 
VISED  diagrams  for  the  cells  and  surface  definitions  were  used  to  generate  input  data 
files  for  the  ITS/ ACCEPT  coupled  electron/photon  Monte  Carlo  program.  The  cosine- 
weighted  current  source  is  a  standard  source  option  in  the  ITS/ACCEPT  pro-am.  It  is 
implemented  by  sampling  the  incident  polar  cosine  of  the  source  particle  as  follows: 

With  fi  defined  as  in  Eq.  1  (above),  then 

M  =  (2) 

where  £  is  a  pseudo-random  number  uniformly  distributed  between  0  and  1.  This  is 
equivalent  to  sampling  /z2  =  £  (3) 

which  then  corresponds  to  the  weighting  of  the  sampled  cosine,  jjl,  by  /j.  itself.  It  can  be 
readily  shown  that  for  the  weighting  function  f{jj)  =  ju  the  cosine-weighted  source  has, 
for  a  negatively  or  inward-directed  normal,  an  average  cosine  given  by 

o 

<M>=-0 - =  (4) 

J  f(/u)dfi 
- 1 

and  an  average  incident  polar  angle  given  by 

o 

Jeos 

<0>==s—a - -J,  (5) 

jf(ju)dfi 


a  necessary  condition  for  uniform  angular  distribution  (isotropic  flux)  with  respect  to  the 
normal  to  the  surface  of  impingence. 

The  cosine-weighted  source  option  is  coded  into  the  original  version  of 
ITS/ACCEPT  for  source  particles  incident  on  a  plane  surface.  To  apply  the  cosine- 
weighted  current  source  methodology  to  the  dome  geometry,  we  modified  the  dome 
source  option  that  we  had  previously  written  for  the  ITS/ ACCEPT  code  [8].  In  Figure  5 
we  define  the  (x,y,z)  axes  to  be  the  inertial  reference  frame  for  the  standard  velocity 
direction  cosines  (a,P,y)  as  defined  in  ITS/ACCEPT,  i.e.  in  terms  of  the  polar  and 
azimuthal  angles  $  and  <p  of  the  particle  trajectory, 

«  =  sin#  cos  (j), 

p  =  sin#  sin  (6) 

y  =  cos#. 
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The  primed  coordinate  system  is  defined  by  a  rotation  through  angle  ©  about  the  x-axis, 
as  shown  in  Figure  5.  The  velocity  direction  cosines  {d,P\'/)  of  the  source  electron  in 
the  primed  system  are  correspondingly 

d  =  sin#' cos  0', 

/?'  =  sin  &  sin  0',  (7) 

/  =  cos  #'. 

The  procedure  to  determine  the  orientation,  expressed  in  the  unprimed  or  inertial  system, 
of  the  individual  source  electrons  is  as  follows: 

1)  choose  the  point  of  particle  incidence  (x0,yo,z0)  to  specify  the  rotation  angle  0  by 
random  uniform  sampling  on  the  hemispherical  surface  (Figure  5); 

2)  sample  cos#'  =  ,  where  as  before,  £  is  a  pseudo-random  number  uniformly 

distributed  between  0  and  1; 

3)  sample  0'  uniformly  on  the  interval  (0, 2 n) ; 

4)  compute  the  direction  cosines  (a,p,y)  as 


1 

0 

0 


0 


0 


ivi 


cos© 

-sin© 


sin© 

cos© 


P 

/ 


(8) 


Figure  5.  Angular  orientation  diagram  for  cosine- weighted  source  incident  on 
CRRES  dosimeter  dome  at  (x0ly0,z0J.  0  is  the  angle  of  rotation  about  the  x-axis. 

#'is  the  polar  angle  between  the  source  electron  direction  and  the  z  ’  axis. 

ITS/ACCEPT  was  modified  (subroutine  HIST)  to  perform  the  above  source  sampling 
procedure.  The  code  modifications  are  shown  in  Appendix  1 . 

Two  sets  of  17  ITS/ACCEPT  runs,  107  histories  each,  were  made  for  the  17 
source  energies  (see  Table  1,  column  1)  listed  in  [1].  Comparisons  of  our  results  with  [1] 
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are  shown  in  Table  1 .  In  our  ITS/ACCEPT  and  MCNPX  runs,  the  calculations  of  energy 
deposition  were  made  using  two  methods;  1)  the  standard  method  of  ACCEPT  and 
MCNPX  -  (energy  entering  cell  -  energy  exiting  cell)[Table  1,  columns  2,4,7];  and  2) 

summation  over  several  energy  bins,  #(£,.)£,.  where  N(Ef)  is  the  pulse  height  in 

I 

energy  bin  i  and  Et  is  the  average  energy  of  the  bin  [Table  1,  columns  3,5,6, 8],  This 
second  method  was  employed  in  [1],  because  the  direct  energy  deposition  calculation  had 
apparently  not  been  installed  in  the  MCNP  version  available  at  the  time  (1993). 

The  first  set  of  runs  employed  the  source  sampling  procedure  described  above.  A 
sample  input  file  for  ITS/ ACCEPT  is  listed  in  Appendix  2.  For  comparison  purposes,  the 
runs  in  the  second  set  were  made  for  the  same  energies,  but  without  cosine  source 
weighting.  The  data  shown  in  column  6  of  Table  1  are  the  results  given  in  [1]  adjusted  to 
correspond  to  the  energy  deposition  units  (MeV)  produced  by  the  ITS/ACCEPT  and 
MCNPX  calculations.  In  [1],  the  energy  deposition  results  are  given  in  units  of  MeV- 
cm-steradian  per  source  electron,  with  the  solid  angle  Q  factor  for  the  CRRESJDl 
geometry  taken  as  34.39  cm2-steradian.  The  data  of  column  6  are  the  raw  results  given  in 
[1]  divided  by  this  solid  angle  factor. 

In  Figure  6  we  plot  the  silicon  dosimeter  (cell  1,  Figure  3)  flux  response  vs. 
source  electron  energy  and  dose  response  vs.  source  electron  energy  curves  obtained  with 
our  calculations  on  the  same  graphs  as  the  Ref.  1  data  for  comparison.  The  flux  and  dose 
response  functions  are  defined  as 

Flux  Response  =  Q]T  N(Ei ) ,  Dose  Response  =  Q]T  EiN(Ej ) , 

*  I 

with  their  respective  units  of  cm2-steradian  and  MeV-cm2-steradian  per  source  electron. 
The  designations  LOLET  and  HILET  refer  to  the  contributions  to  the  response  functions 
from  electrons  with  energy  below  and  above  1  MeV  in  the  dosimeter  (cell  1,  Figures), 
respectively. 

Table  2  shows  a  comparison  of  ITS/ACCEPT  and  MCNPX  energy  deposition 
results  in  all  cells  (Figure  3)  for  the  four  source  energies  (2,  5,  10,  25  MeV)  for  which 
both  programs  were  run. 
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Electron  Energy  (MeV) 


0  5  10  15  20  25 


Electron  Energy  (MeV) 

RFF  1  *  Anrhamnnnnh  A  finvtnn 

Figure  6.  Comparison  of  Flux  and  Dose  Response  functions  calculated 
with  ITS/ACCEPT[5]  and  as  reported  in  [1]. 


10 


o 


oo 


> 

m 

a 

tn 

<N 

11 

> 

a 

& 

MCNPX 

2.729E-04 

9.013E-04 

1.862E+00 

2.360E-01 

7.020E-02 

1.029E-01 

5.948E-02 

1.209E-01 

4.130E+00 

4.116E-02 

1.622E-01 

4.866E-01 

7.273  E+00 

m 

& 

m 

m 

u 

s£ 

D 

O 

00 

ACCEPT 

2.215E-04 

7.458E-04 

2.205E+00 

1.732E-01 

6.147E-02 

7.540E-02 

1.342E-02 

1.204E-01 

O 

I 

SO 

3.353E-02 

1.216E-01 

3.895E-01 

7.355E+00 

SOURCE  ENERGY  =  10  MEV 

MCNPX 

2.828E-04 

9.973E-04 

1.836E+00 

2.447E-01 

7.610E-02 

1.057E-01 

3.348E-02 

<N 

O 

m 

s 

o 

1.397E+00 

4.343E-02 

1.662E-01 

5.142E-01 

4.488E+00 

H 

CU 

m 

u 

u 

< 

2.380E-04 

O 

i 

m 

t — 

On 

Os 

2.124E+00 

1.943E-01 

6.787E-02 

8.759E-02 

1.665E-02 

6.323E-02 

1.530E+00 

3.963E-02 

1.450E-01 

4.907E-01 

4.760E+00 

SOURCE  ENERGY  =  5  MEV 

MCNPX 

2.898E-04 

1.033E-03 

1.687E+00 

o 

I 

m 

cn 

<N 

7.564E-02 

9.467E-02 

1.738E-02 

3.233E-02 

3.01  IE-01 

3.726E-02 

1.457E-01 

4.224E-01 

3.046E+00 

ACCEPT 

2.668E-04 

9.164E-04 

1.875E+00 

1.908E-01 

8.126E-02 

8.851E-02 

2.737E-03 

1.559E-02 

2.902E-01 

4.812E-02 

1.440E-01 

© 

1 

m 

00 

00+398  IT 

Source  ENERGY  =  2  MEV 

MCNPX 

1.790E-04 

6.471E-04 

1.225E+00 

1.054E-01 

3.273E-02 

2.445E-02 

9.523 E-03 

5.399E-03 

5.721E-03 

s 

<N 

CN 

<N 

O 

m 

so 

CN 

CN 

2. 1 26E-02 

1.455E+00 

ACCEPT 

1.975E-04 

6.491E-04 

00+380TI 

6.332E-02 

4.008E-02 

1.204E-02 

6.120E-06 

4.050E-05 

3.602E-03 

2.292E-03 

7. 183  E-03 

1.394E-03 

1.439E+00 

MCNPX 
CELL# 
(FIGURE  3) 

SO 

+ 

so 

O 

*•— i 

so 

<N 

"3* 

t 

TOTAL 

11 


2.2  HEP-  Electron  Transport  Modeling 

Extensive  testing  was  performed  on  the  ITS/ ACCEPT  input  file  written[9]  for  the 
HEP  instrument.  This  file  contains  a  complete,  blueprint  matching,  geometry  description 
for  the  in-flight  model.  It  was  felt  that,  due  to  the  complicated  geometry  structure  of  the 
model,  it  would  be  advisable  to  run  an  extensive  series  of  Monte  Carlo  calculations  for 
model  validation.  To  that  end  four  sets,  each  consisting  of  ten  source  positions,  of 
(50000  history)  electron  transport  runs  were  made.  The  first  of  these  consisted  of  pencil 
beam,  15  MeV,  sources  normally  incident  on  the  “top”  (x  =  5.0  cm)  at  ten  uniformly 
distributed  positions  along  z.  The  second  and  third  sets  were  similarly  distributed  along  z 
but  were  positioned  on  the  side  (y  =  4.0  cm)  and  bottom  (x  =  -3.0  cm).  The  fourth  set 
consisted  of  calculations  made  with  a  6-component  energy  (max.  energy  =  20  MeV) 
spectrum,  point  isotropic  source  embedded  at  ten  uniformly  spaced  positions  along  the 
central  (z)  axis.  In  all  cases  run,  there  were  no  “lost  particles”  attributable  to  geometry 
specification  errors. 

3.  PROTON  TRANSPORT 

3.1  CEASE  -  MCNPX  Proton  Transport  Modeling 

Two  proton  transport  runs  were  made  with  the  CEASE  telescope  model  [10]  for 
MCNPX,  as  requested  by  the  Air  Force  sponsor,  for  9  MeV  proton  flat  disk  sources, 
located  at  z  =  0.457  cm  on  the  telescope  axis,  at  40  degree  slant  incidence  with  respect  to 
the  z-axis.  Partial  listings  of  the  run  outputs  are  shown  in  Appendix  3.  The  geometry 
portions  of  these  files,  previously  listed  in  [1 1]  were  omitted  for  brevity. 

3.2  Grazing  Angle  Proton  Scattering  Calculations 

An  extensive  series  of  proton  transport  calculations  were  made  using  MCNPX  for 
protons  beams  incident  on  aluminum  and  iridium  (6x5x1  cm)  at  shallow  (grazing) 
angles  of  incidence  ( 9in  =90.1, 90.5, 91.0°,  <|>  =  90.0°  as  shown  in  Figure  7). 


Figure  7.  Problem  geometry  -  proton  beam  incident  on  material  slab 

The  motivation  for  this  work  was  an  investigation  into  the  occurrence  of  unexpected 
radiation  damage  to  one  of  the  Chandra  X-Ray  Telescope  cameras  by  grazing  incidence 
protons.  Standard  Monte  Carlo  codes,  such  as  Geant4[12],  the  program  used  to  predict 
radiation  damage  in  the  Chandra,  produced  erroneous  results.  In  our  calculations  we  set 
out  to  investigate  the  treatment  by  standard  Monte  Carlo  programs  such  as  MCNPX  of 
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the  energy  loss  and  relationship  between  the  beam  angle  of  incidence  and  die  emergent 
angular  distribution  (from  the  top  surface  or  vacuum  interface  of  a  material  medium)  for 
protons  with  grazing  incident  angles.  To  accomplish  this,  we  ran  MCNPX  for  0.25  MeV 
protons  incident  on  aluminum  and  0.05,  0.25,  0.5,  1.0  MeV  protons  incident  on  iridium 
slabs  at  grazing  angles.  All  of  the  Monte  Carlo  calculations  were  made  using  10®  case 
histories  to  ensure  statistical  reliability.  High  angular  resolution  emergent  proton  spectra 
were  obtained,  and  plots  made,  for  energy-angle,  total  energy,  and  azimuthal 
distributions. 


The  MCNPX  proton  transport  results  were  analyzed  from  several  aspects. 
Programs  were  written  to  convert  the  raw  Monte  Carlo  data  to  histograms  detailing: 

„  /to -so-)’  - 


(1)  probability  density  Pt 


E„ 


of  emergent  protons  with 


exit  polar  angle  0OU, ,  for  every  fixed  incident  angle  0in ,  for  every  source  energy 
Ein  and  exit  energy  bin  width  A£OT,  =  0.05£fe .  Representative  plots  of  the 
probability  density  function  for  0.25  MeV(=£/„)  protons  incident  on  iridium  with 
incident  angle  90.5  degrees(=  6in ),  are  shown  in  Figure  8.  The  five  plots 

represent  the  emergent  proton  probability  densities  for  the  five  emergent  energy 
intervals  0-12.5,  12.5-25.0, 25.0-37.5, 37.5  50.0,  50.0-62.5  keV. 


(2)  emergent  proton  pulse-height  distribution  histogram  plots  N(Eout,0out\Ein,0jn)\ 
Representative  plots  are  shown  in  Figures  9  and  10  for  0.25  MeV(=£jn)  protons 
incident  on  iridium  with  incident  angle  91.0  degrees(=  0in ).  The  emergent  proton 
pulse-height  distributions  were  calculated  using  5  degree  coarse  resolution 
emergent  angle(0ou()  bins  in  the  range  O<0gut  <85°  (Figure  9)  and  0.1  degree 

high  resolution  bins  in  the  range  85°  <  0Bu!  <  90°  (Figure  10).  The  pulse-height 

distributions,  expressed  as  the  number  of  emergent  protons  per  energy-angle  bin, 
per  incident  proton,  are  normalized  to  unit  proton  incidence. 

(3)  emergent  proton  pulse-height  energy  spectra  N(Eout;Ejn,0in) 
C—  ^  N(Eout ,  0OU, ,  Ejn ,  0in )  ] .  The  plot  shown  in  Figure  11  is  the  pulse-height 

®out 

energy  spectrum  of  emergent  protons  resulting  from  a  0.25  MeV  proton  beam 
incident  on  iridium  with  angle  91 .0  degrees(=  0in ). 

(4)  emergent  proton  azimuthal  ($,„, )  distribution.  The  plot  shown  in  Figure  12  is  the 
distribution  of  (f>  out  for  0.5  MeV  protons  incident  on  iridium  with  0in  —  91.0 
degrees.  The  emergent  azimuthal  distribution  is  strongly  peaked  in  the  incident 
azimuthal  direction,  $„= 90.0  degrees.  This  was  found  to  be  true  for  all  15  cases 
studied. 
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Additionally,  plots  were  made  for  each  incident  energy-angle  combination,  of  the 
maximum  value  of  the  probability  density  Pmsx{Eout',Ein,6in)  (Figure  13)  and  its 


corresponding 


ratio  vs.  Eout  (Figure  14). 


The  results  obtained  with  the  MCNPX  model  do  not  exhibit  specular,  or  near- 
specular,  reflection[13]  as  may  have  been  expected  from  experimental  results  [14,15]. 
Preliminary  conclusions  also  reveal  that  the  computed  proton  energy  loss  resulting  from 
scatter  of  grazing  incidence  protons  is  an  overestimation  and  are  significantly  larger  than 
experimental  observation  [14,15].  In  MCNPX  and  Geant4,  proton  energy  loss 
calculations  in  the  energy  range  studied  are  based  on  the  assumption  that  the  continuous 
slowing-down  approximation  applied  in  conjunction  with  the  bulk  material  stopping 
power  is  appropriate.  Pfandzelter  et  al.  [15]  have  derived,  from  a  study  of  the  physics  of 
surface  potentials,  expressions  for  the  effective  stopping  power  that  applies  in  this 
situation. 

Sample  MCNPX  output  is  shown  in  Appendix  4  for  a  250  keV  proton  beam 
incident  on  iridium  at  0.1°  (0iH  =  90.1°  as  shown  in  Figure  7).  The  energy  and  angular 

bin  structures  of  output  emergent  proton  current  tallies  shown  are  as  defined  in  items  1-3, 
above.  For  these  investigations,  a  specialized  source  subroutine,  listed  in  Appendix5, 
was  written  for  MCNPX.  The  azimuthal  {(pout  )  distribution  of  the  emergent  proton 
current  was  extracted  from  PTRAC  files  produced  in  the  MCNPX  runs.  A  program 
(PHILOOK)  was  written  to  extract  the  azimuthal  distribution  data  and  convert  them  to 
the  histogram  form  exemplified  by  the  polar  plot  of  Figure  12. 
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Figure  9.  Emergent  proton  pulse-height  d: 
protons  incident  on  Iridium  slab  at  0in  =9 
increments. 


Figure  1 0.  Emergent  proton  pulse-height  distribution  N(E0Ut ,  $oul ;  Ein ,  6in )  for  250  keV 
protons  incident  on  Wdium  slab  at  6in  =  91.0° .  Results  are  shown  for  85°  <  60ut  <  90°  in 
0.1°  increments. 


0.000  0.025  0.050  0.075  0.100  0.125  0.150  0.175  0.200  0.225  0.250 

Eoui MeV) 


Figure  1 1 .  Pulse-height  energy  spectrum  N(E0Ut ;  Ein ,  6in )  for  protons  emerging  from  iridium 
slab.  E in-250  keV,  0,n=9O.5o. 
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Eoul  MeV) 

Figure  1 3 .  Maximum  value  of  the  probability  density  (Eoul ;  Ein ,  0^ )  vs.E0Ut  for  protons 


of  energy  Ein=l. 0  MeV  incident  on  iridium  at  6in  =  90.1°. 


4.  SUMMARY 


During  the  period  covered  by  this  report,  the  technical  activity  and  progress 
achieved  consisted  primarily  of:  1)  modeling  and  verification  of  earlier  calculations  of 
electron  transport  and  energy  deposition  in  CRRES  dosimeters;  2)  model  validation  of 
electron  transport  simulations  in  the  HEP  instrument;  3)  modeling  proton  transport  in  the 
CEASE  instrument;  4)  the  beginning  of  a  study  of  the  scattering  of  grazing  incidence 
protons  from  material  surfaces;  and  5)  providing  computer  code  enhancements  and 
advice  on  Monte  Carlo  simulation  code  implementation  to  AFRL, 


ft 
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Appendix  l 


SUBROUTINE  HIST  IS  CALLED  BY 


SUBROUTINE  HIST  CALLS 
INTRINSIC  FUNCTIONS 


EXTERNAL  FUNCTIONS 


SQRT, 

REAL 


ITS 


RANF 


CLASS,  ECHOS,  EHIST,  TIMER, 

RANINT,  RANSAV 

ZONE 

FOLD,  ZONEA 
PLTDAT 


[CYLTRAN) 

PHI  ST 


Subroutine  HIST  for  Cosine-Weighted  Dome  Source 

SUBROUTINE  HIST  HIST 

***********************************  *  *  +  *  *  ****  *  *  *  *  *  *  *  *  *  *  *  *  **  *  ^  ^  *  **  ^  *|jjgij! 

HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 
HIST 

************************************************** *******************+HIST 


(CYLTRAN) 

(ACCEPT) 

(M-CODES) 


ORIGINATION  DATE 
LAST  MODIFIED 


16  JAN  68. 
30  MAY  91 


FUNCTION 

THIS  PROGRAM  SAMPLES  PHASE  SPACE  PARAMETERS  FOR 
SOURCE  PARTICLES.  SUBSEQUENTLY  CALLS  EITHER  EHIST  OR 
PHI ST.  RETRIEVES  "BANKED"  ELECTRONS  AND  CALLS  EHIST. 
TALLIES  PULSE  HEIGHT  DISTRIBUTION. 


***  COMMON  BLOCKS  CNSTNT,  PARAMS,  OUT,  CALC,  XPED,  STOR,  STTS 
(PAREM) -ACCEPT 

C$  LIST (S=0) 

CDIR$  NOLIST 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

SAVE 


HIST 

HIST 

HIST 

HIST 

CNSTNT 

CNSTNT 


00007 

00009 

00010 

00011 

00012 

00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 

00021 

00022 

00023 

00024 

00025 

00026 

00027 

00028 

00029 

00030 

00031 

00032 

00033 

00034 

00035 

00036 

00037 

00081 

00082 


C 


PARAMETER  (CT1M12  =  1.0D-12,  CT1EM8  =  1.0D-8,  CT1EM7 


CNSTNT  00083 
1.0D-7)  CNSTNT  00084 


CNSTNT  common  block  is  identical  to  that  given  in  Reference  3 
-  listing  omitted  for  brevity. 


PARAMETER  (CCOHIM=57 . 031547D0 ,  CCOHMX=80 . 654788D0) 

C - - - - — - - - - - - 


CNSTNT  00140 
PARAMS  00002 
PARAMS  00003 


PARAMS  common  block  is  identical  to  that  given  in  Appendix 
2,  Reference  4. 


c 


PARAMETER  (NLAST  =  50) 

PARAMETER  {  INUMK  =  3,  INGP  =  INMT) 

COMMON  /EXTSORC/  IRECTS ,  IDISKS ,  XLOWS ,  XHIGHS ,  YLOWS ,  YHIGHS ,  ZLOWS 
$  ZHIGHS ,  XCENT ,  YCENT ,  ZCENT ,  XCIR ,  YCIR ,  2CIR,  KPERPYZ ,  KPREPXZ ,  KPERPXY , 
$  IDOME ,  RDOME 
LOGICAL  RRKILL,  FLMTEL 
COMMON  /OUT/ 


PARAMS  00086 
PARAMS  00089 
PARAMS  00103 


OUT  00002 

OUT  00003 


25 


o  o  o  on 


1  FLMTEL ( INGP ) 


OUT  common  block  is  identical  to  that  given  in  Reference  3  - 
listing  omitted  for  brevity. 


$  FLESCP,  FLNOK,  FLNEL,  FLBAD,  FLINIT,  FLGSEC  OUT  00095 

CALC  00002 

COMMON  /CALC/  CALC  00003 

2  CALC  common  block  is  identical  to  that  given  in  Reference  3  -  34 


listing  omitted  for  brevity. 


COMMON  /XPED/ 

1  DETOUR (INMT) ,  RHO ( INMT ) ,  MT,  MTP ,  MTP0 

LOGICAL  DMPFLG ,  FLMC 
DOUBLE  PRECISION  IRSAV 

COMMON  /STTS/  IB,  NB,  NSORS,  IBT,  BC 
$  ,  IHIST,  IRSAV,  KPUT,  FLMC 


CHARACTER *3  OTYPE(IO),  OBODY 
LOGICAL  FLDBG ,  FLDBGL 
COMMON  /  PAREM/ 

$  XB (3 )  ,  WT ( 3 )  ,  RIN,  ROUT, 

$  FLDBG,  IRPRIM,  ICALL,  LSURF, 

$  KLOOP,  LOOP,  ITYPE,  FLDBGL 

COMMON  /PAREMO/  OTYPE 


$  IHSTRY, COFSRC (6,25,10), COFLEGY (6,25), COFLEG ( 6 ) 

C$  LIST ( S=1 ) 

CDIR$  LIST 

COMMON  /STOR/ 


1  CTHS(NLAST), 

2  LBS (NLAST)  , 

$  ,XS (NLAST), 

3  CPHS (NLAST), 

4  ,LBCS (NLAST) 

EXTERNAL  RAN 


TS (NLAST) , 
NTS (NLAST) 
YS (NLAST) , 
SPHS (NLAST) 


WS (NLAST) ,  ZS (NLAST) ,  IPRS (NLAST) , 
STHS (NLAST) , 


CIMAX  =  IMAX 
IF  (FLSPEC)  THEN 
TAV  =  CZERO 
ELSE 

TAV  =  CIMAX*TIN 
END  IF 

CALL  RANINT  ( IRA) 

IF  (IB  .EQ.  1)  INRAN  =  IRA 
DO  130  1=1,  IMAX 
IHSTRY=I 

DO  1301  JJJ=1,10 
EDPR ( JJJ) =0 . 

EDNK ( JJJ) =0 . 

EDSC (JJJ) =0 . 

EDTL (JJJ) =0 . 

IHIST  =  I 

MODTMJ  =  MIN  ( 10000 , IMAX) 

IF ( I . EQ . MODTMJ* ( I /MODTMJ ) )  THEN 


uuxi9 

XPED 

00002 

XPED 

00003 

XPED 

00010 

XPED 

00012 

STTS 

00002 

STTS 

00010 

STTS 

00017 

STTS 

00018 

STTS 

00019 

PAREM 

00002 

PAREM 

00003 

PAREM 

00004 

PAREM 

00008 

PAREM 

00009 

PAREM 

00013 

PAREM 

00014 

PAREM 

00015 

PAREM 

V, 

00016 

HIST 

00047 

HIST 

00048 

STOR 

00002 

STOR 

00003 

STOR 

00004 

STOR 

00006 

STOR 

00007 

STOR 

00009 

HIST 

00050 

RANNUM 

00003 

HIST 

00089 

HIST 

00090 

HIST 

00091 

HIST 

00092 

HIST 

00093 

HIST 

00094 

HIST 

00095 

HIST 

00096 

HIST 

00097 

HIST 

00098 

HIST 

00101 

HIST 

00103 
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non  o  o  ooooo 


CALL  TOTTIM  ( XTM J ) 

WRITE  { * ,  1  ( / 1  '  HISTORY'  ',18,  1  %  ELAPSED  MINUTES  ”,  F10 . 2 )' ) 
lI,XTMJ/60. 

ENDIF 

W  =  CONE 
CWCF  =  W 
LAST  =  0 


CALL  RANSAV  { IRSAV) 


. . .  SOURCE  ENERGY 


14 

16 


20 


IF  (FLSPEC)  THEN 
RA  =  RAN (IRAN) 

DO  14  JHIST  =  2, JSPEC 

IF  (  RA.  GT.  SPECIN (JHIST)  )  GO  TO  16 
CONTINUE 

T  =  ESP(JHIST-l)  -i-  (  RA  -SPECIN ( JHIST- 1)  )*(  ESP(JHIST) 
-  ESP  ( JHIST-1)  )/(  SPECIN  (JHIST)  -  SPECIN  ( JHIST- 1)  ) 
TAV  =  TAV  +  T 

IF  (  (FLESRC  .AND.  (T  .GT.  TOUT  ))  .OR. 

(.NOT.  FLESRC  .AND.  (T  .GT.  TPCUT)  )  )  THEN 
GO  TO  20 

ELSE 

NTREJ  =  NTREJ  +  1 
TREJ  =  TREJ  +  W*T 
GO  TO  1299 
END  IF 
END  IF 
T  =  TIN 
NT  =  NTFST 


CALL  CLASS  (T,NT) 

I F ( I DOME . EQ . 0 ) THEN 

SOURCE  DIRECTION 

IF  (ICTH  .EQ.  2)  THEN 
RA  »  RAN  (IRAN) 

COM  =  CTHIN+  RA*  (CONE-CTHIN) 

ELSE  IF  (ICTH  .EQ.  3)  THEN 
RA  =  RAN(IRAN) 

COM  =  SQRT  (CTHIN+RA*  (CONE-CTHIN)  ) 


LAHEY 

LAHEY 

LAHEY 

LAHEY 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 

HIST 


C 

C 


ELSE  IF  (ICTH  .EQ.  1)  THEN 

HIST 

CTH ( 1 )  =  CTSR 

HIST 

STH ( 1 )  =  STSR 

HIST 

CPH{1)  =  CPSR 

HIST 

SPH(l)  =  SPSR 

HIST 

GO  TO  69 

HIST 

END  IF 

HIST 

HIST 

IF  (CTSR  .EQ.  CONE)  THEN 

HIST 

CTH ( 1)  =  COM 

HIST 

STH (1)  =  SQRT (CONE -COM* COM) 

HIST 

RA  =  RAN(IRAN) 

HIST 

JAZ  =  RA*C36Q 

HIST 

CPHU)  =  CCH  ( JA2+1) 

HIST 

SPH(l)  *  SCH(JAZ+1) 

HIST 

ELSE 

HIST 

HIST 

CALL  FOLD ( CTSR ,  STSR , CPSR ,  SPSR, COM, CTH { 1 )  ,  STH { 1 )  ,CFH(1)  ,SPH(1)  JHIST 

— - - - — - — - - - - - - - — HIST 

END  IF  HIST 

END  IF 

_  HIST 


J 


00019 

00020 

00021 

00022 

00105 

00106 

00107 

00108 

00109 

00110 

00111 

00112 

00113 

00114 

00115 

00116 

00117 

00118 

00119 

00120 

00121 

00122 

00123 

00124 

00125 

00126 

00127 

00128 

00129 

00130 

00131 

00132 

00133 

00134 

00135 

00136 

00137 

00138 

00139 

00140 

00141 

00142 

00143 

00144 

00145 

00146 

00147 

00149 

00150 

00151 

00153 

00154 

00155 

00156 

00157 

00159 

00160 

00161 

00162 

00163 

00165 

00172 

00173 

00174 

00176 

00177 
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o  o 


69 


SOURCE  POSITION 

HIST 

00178 

HIST 

00179 

IF  (SORCIN  .NE.  CZERO)  THEN 

HIST 

00198 

RA  =  RAN  (IRAN) 

HIST 

00199 

R  =  SQRT (RA) * SORCIN 

HIST 

00200 

RA  =  RAN  ( IRAN) 

HIST 

00201 

JAZ  =  RA*C360 

HIST 

00202 

SCHR  =  SCH ( JAZ+1 ) *R 

HIST 

00203 

CCHR  =  CCH (JAZ+1 ) *R 

HIST 

00204 

IF  (IDISKS  .EQ.  0)  THEN 

X  =  XSR  +  CCHR  * W 1 X + S  CHR  *  W2  X 
Y  =  YSR+CCHR*W1Y+SCHR*W2Y 

Z  =  ZSR+CCHR*WlZ+SCHR*W2Z 

ELSE 

I F ( KPERPXY . EQ . 1 ) THEN 
X  =  XCENT  +  CCHR 

Y  =  YCENT  +  SCHR 
Z  =  ZCENT 

END  IF 

IF  ( KPERPXZ . EQ . 1 ) THEN 
X  =  XCENT  +  CCHR 

Y  =  YCENT 

Z  =  ZCENT  +  SCHR 
END  IF 

IF (KPERPYZ . EQ . 1 ) THEN 

X  =  XCENT 

Y  =  YCENT  +  CCHR 

Z  =  ZCENT  +  SCHR 


END  IF 
END  IF 
ELSE 

HIST  00208 

HIST  00209 

HIST  00210 

HIST  00211 


IF  (IRECTS  .EQ.  0  .AND.  IDOME . EQ . 0 ) THEN 
X  =  XSR 
Y  =  YSR 
Z  =  ZSR 
ELSE 

IF ( IRECTS . NE . 0 ) THEN 


RRAA1  =  RAN  ( IRAN) 


RRAA2  =  RAN  {I RAN) 

IF  (KPERPXY  .EQ.  1)  THEN 

X  =  XLOWS  +  RRAA1* (XHIGHS-XLOWS ) 
Y  =  YLOWS  +  RRAA2* ( YHIGHS-YLOWS ) 
Z  =  ZLOWS 


END  IF 

IF  (KPERPXZ.  EQ.  1)  THEN 

X  =  XLOWS  +  RRAA1* (XHIGHS-XLOWS) 
Y  =  YLOWS 

Z  =  ZLOWS  +  RRAA2MZHIGHS- ZLOWS) 
END  IF 

IF  (KPERPYZ  . EQ . 1 )  THEN 
X  =  XLOWS 


C 
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no  o  o  non  o  o  ono  norm 


Y  =  YLOWS  +  RRAAl* (YHXGHS-YLOWS) 

Z  =  ZLOWS  +  RRAA2* (ZHIGHS- SLOWS) 

END  IF 

^ — —— i  i  — — ^ — — — — — — 

C  DOME  SOURCE 

IF { IDOME . HE . 0 ) THEN 
STHDM=RAN  { IRAN } 

CTHDM=SQRT ( 1 . - STHDM *  STHDM ) 

PPHDM=C2PI *RAN ( IRAN) 

CPPHDM=COS (PPHDM) 

SPPHDM=SIN { PPHDM) 

ALDM = STHDM  *  CP  PHDM 
BTDM” STHDM  * SPPHDM 
X=RDOME*ALDM 
Y = RDOME  *  BTDM 
Z^RDOME  *CTHDM 

106  PHSDM=C2PI*RAN  (IRAN) 

C  COSINE -WEIGHTED  SOURCE 

CTHSDM=  -SORT  {RAN  { IRAN}  } 

STHSDM=SQRT { 1 , -CTHSDM* CTHSDM ) 

S PHSDM- SIN { PHSDM ) 

CPHSDM=COS (PHSDM) 

UUUD=STHSDM*CPHSDM 
WVD=STH3DM*  SPHSDM 
WWWD=CTHSDM 

AAAD=SQRT(UUUD**2+WVD**2*WWWD**2) 
UUUD=UUUB/AAAD 
VWD=VWD/AAAD 
WWWD=WWWD  /  AAAD 
ALPDD=UUUD 

BETDD=WVD*CTHDM+WWWD*  STHDM 
GAMDD=  -WVD*  STHDM +WWWD  *  CTHDM 

PRODD- ALDM  *  ALFDD-f  BTDM  *BETDD+  CTHDM*GAMDD 
IF (PROBD.GT.O.OJGO  TO  106 
STH ( 1 } =SQRT ( 1 , -GAMDB**2 ) 

CTH ( 1 } =GAMDD 

IF ( STH ( 1 }  , LT . 1 . E-6 } GO  TO  106 
SPH (1 ) =BETDD/ STH { 1 ) 

CPH { 1 } =ALPDD/STH { 1 ) 

.  END  IF 

BHD  IF - 


END  IF 
C 

XB(1)  =  X 
XB{2)  =  Y 
XB (3)  =  Z 

WT(1)  =  STH { 1 ) *CFH { 1 ) 
WT(2}  =  STH (1) *SPH{1) 
WT(3)  =  CTH ( 1 } 


CALL  ZONE A 

C  - — - - - 

LB  =  IR 
LBCZ  =  IRPRIM 
IPR  =  1 


.  . .  CALL  TRACKING  ROUTINES 


70  IF  {  FLESRC  .OR.  (IPR  .NE.  1}  }  THEN 


PARTICLE  TO  BE  TRACKED  IS  AN  ELECTRON 


IF  {MT  .NE.  MAT (LB) )  THEN 
MT  =  MAT (LB) 

END  IF 


CALL  EH I ST 


ELSE 


PARTICLE  TO  BE  TRACKED  IS  A  PHOTON 


LPCZ  =  LBCZ 


CALL  PHIST(X,Y#Z,LB/CTH{1) ,3TH{1) #CPH{1) ,SPH(1) #T#W# 1) 


END  IF 


HIST 

00212 

HIST 

00213 

HIST 

00220 

HIST 

00221 

HIST 

00222 

HIST 

00223 

HIST 

00224 

HIST 

00225 

HIST 

00226 

HIST 

00227 

HIST 

00228 

HIST 

00229 

HIST 

00230 

HIST 

00232 

HIST 

00233 

HIST 

00234 

HIST 

00235 

HIST 

00236 

HIST 

00237 

HIST 

00238 

HIST 

00239 

HIST 

00240 

HIST 

00241 

HIST 

00242 

HIST 

00248 

HIST 

00249 

HIST 

00250 

HIST 

00251 

HIST 

00252 

HIST 

00253 

HIST 

00254 

HIST 

00255 

HIST 

00262 

HIST 

00265 

HIST 

00266 

HIST 

00267 

HIST 

00269 

HIST 

00270 

HIST 

00271 
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onooo  o  o  no 


.  .  .  REMOVE  SECONDARY  ELECTRONS  FROM  STORAGE  FOR  TRANSPORT 


IF  (LAST  .NE.  0)  THEN 
LB  =  LBS (LAST) 

Z  =  ZS (LAST) 

T  =  TS (LAST) 

NT  =  NTS (LAST) 

CTH(l)  =  CTHS (LAST) 

W  =  WS ( LAST ) 

I PR  =  I PRS ( LAST ) 

X  =  XS ( LAST ) 

Y  =  YS (LAST) 

STH(l)  =  STHS (LAST) 

CPH(l)  =  CPHS (LAST) 

SPH(l)  =  SPHS (LAST) 

LBCZ  =  LBCS (LAST) 

KLOOP  =  KLOOP+1 
LAST  =  LAST-1 
GO  TO  70 
END  IF 

IF  (.NOT.  FLPHD )  GO  TO  1299 


.  . .  SCORE  PULSE-HEIGHT  DISTRIBUTION 


EABST  =  CZERO 

DO  100  LS=LPHDB , LPHDE 

EABST  =  EABST+PHDD ( LS ) 

100  PHDD(LS)  =  CZERO 
DO  110  JS=1,JSMAX 

IF (SMARK ( JS)  .LE.  EABST)  GO  TO  120 
110  CONTINUE 

NPHD  =  NPHD+1 
GO  TO  1299 

120  ABE(JS)  =  ABE ( JS) +CWCF 
1299  IF (NINDV. EQ. 0 ) GO  TO  130 

DO  1298  NIND=1# NINDV 

EDTL  (NIND)  =EDPR  (NIND)  +EDNK  (NIND)  +EDSC  (NIND) 

1298  CONTINUE 

WRITE  (44  )  (EDPR(NIND)  ,  EDNK(NIND)  ,EDSC(NIND)  y  EDTL  (NIND)  ,  NIND 
$  =1, NINDV) 

130  CONTINUE 

C  _ 

CALL  RANSAV ( IRC) 


RETURN 

END 


HIST  00272 
HIST  00273 
HIST  00274 
HIST  00275 
HIST  00276 
HIST  00277 
HIST  00278 
HIST  00279 
HIST  00280 
HIST  00281 
HIST  00283 
HIST  00284 
HIST  00285 
HIST  00286 
HIST  00287 
HIST  00288 
HIST  00289 
HIST  00291 
HIST  00292 
HIST  00294 
HIST  00295 
HIST  00296 
HIST  00297 
HIST  00298 
HIST  00299 
HIST  00300 
HIST  00301 
HIST  00302 
HIST  00303 
HIST  00304 
HIST  00305 
HIST  00306 
HIST  00307 
HIST  00308 
HIST  00309 
HIST  00310 
HIST  00311 
HIST  00312 


HIST  00313 
HIST  00314 

HIST  00315 

HIST  00316 

HIST  00317 

HIST  00318 
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Appendix  2 

CRRES_D1  Sample  Input  File  for  ACCEPT  Run 


TITLE 

15  MEV  DOME  SOURCE  on  CRRES  DOME  1 


**********************  GEOMETRY  *********************** 
GEOMETRY 
*1 

RCC 

*2 

o 

o 

o 

o 

o 

o 

o 

o 

0.00000 

0.04030 

.05093 

RCC 

*3 

o 

o 

o 

o 

-0.07620 

o 

o 

.00000 

0.07620 

.06350 

RPP 

*4 

-.09525  . 

09525  -.09525 

.09525  .00000  0.04030 

RCC 

o 

o 

o 

o 

-0.07620 

o 

o 

.00000 

0.07620 

0.5080 

RCC 

*6 

o 

o 

o 

o 

-0.07620 

o 

o 

.00000 

0.07620 

1.22936 

RCC 

*7 

o 

o 

o 

o 

-0.07620 

o 

o 

.00000 

0.07620 

1.46812 

RCC 

*8 

o 

o 

o 

o 

-0.10668 

o 

o 

.00000 

0.03048 

0.32385 

RCC 

*9 

o 

o 

o 

o 

-0.10668 

o 

o 

.00000 

0.03048 

1.22936 

RCC 

*10 

o 

o 

o 

o 

-0.10668 

o 

o 

.00000 

0.03048 

1.46812 

RCC 

*11 

0.0 

0.0 

-0.15748 

o 

o 

.00000 

0.05080 

0.32385 

RCC 

*12 

0.0 

0.0 

-0.15748 

o 

o 

.00000 

0.05080 

1.22936 

RCC 

*13 

0.0 

o 

o 

-0.15748 

0.0 

.00000 

0.05080 

1.46812 

RCC 

*14 

o 

o 

0.0 

-0.31496 

0.0 

.00000 

0.15748 

1.22936 

RCC 

*15 

o 

o 

o 

o 

-0.31496 

o 

o 

.00000 

0.15748 

1.46812 

RCC 

*16 

o 

© 

o 

o 

-0.94996 

o 

o 

.00000 

0.635 

1.11125 

RCC 

*17 

o 

o 

o 

o 

-0.94996 

o 

o 

.00000 

0.635 

1.22936 

RCC 

*18 

o 

o 

o 

o 

-0.94996 

o 

o 

.00000 

0.635 

1.46812 

RCC 

*19 

o 

o 

o 

o 

-1.58496 

o 

o 

.00000 

0.635 

1.11125 

RCC 

*20 

o 

o 

o 

o 

-1.58496 

0.0 

.00000 

0.635 

1.22936 

RCC 

*21 

o 

o 

o 

o 

-1.58496 

o 

o 

.00000 

0.635 

1.46812 

RCC 

*22 

o 

o 

o 

o 

o 

o 

o 

o 

.00000 

1.01981 

1.01981 

RCC 

*23 

o 

o 

o 

o 

o 

o 

o 

o 

.00000 

1.22936 

1.22936 

RCC 

*24 

SPH 

*25 

SPH 

*26 

SPH 

*27 

o  o  o  o 

o  o  o  o 

0.0 

0.0 

0.0 

0.0 

0.0  0.0 

0.0  1.01981 

0.0  1.22936 

0.0  1.46812 

.00000 

1.46812 

1.46812 

RCC 

*28 

SPH 

*29 

0.0 

0.0 

o  o 

o  o 

-2.58496  0. 

0.0  1.005 

0  0.0  1 

.0  1.46812 

RCC 

*30 

SPH 

*31 

SPH 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o 

o 

o 

•  in 
o 

o 

o 

o  *  o 

*  o 
o 

o 

O  H 

1.005  1.005 

0 
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END 

*SI 

Z01  +1 
Z02  +3  -1 
♦VOID 

Z03  +29  +28  -1  -3 
*AL 

Z04  +25  -24  +22 
Z05  +5  -4  -2 
♦VOID 

z06  +6  -5  -4  -2 
♦AL203 

Z07  +4  -2 
♦VOID 
Z08  +2 
*AL 

Z09  +8  -7 
♦VOID 

zlO  +9  -8  -7 
Zll  +7 

zl2  +12  -11  -10 
*AL 

Z13  +11  -10 
*NI 

Z14  +10 
*AL 

Z15  +13 
♦VOID 

zl6  +14  -13 
zl7  +17  -16 
♦AL 

zl8  +16  -15 
zl9  +19  -18 
*W-Fe 
Z20  +15 
z21  +18 
♦AL 

z22  +20  -19 
♦VOID 

Z23  +24  -28  +21 
Z24  +23  -25  +26 
z25  +23  -26 
Z26  +27 

Z27  +29  -1  -2  -3  -4  -5  -6  -7  -8  -9  -10  -11  -12  -13  -14  -15 
-16  -17  -18  -19  -20  -21  -22  -23  -24  -25  -26  -27  -28 
Z28  +30  -29 
END 

♦MATERIAL 

1 

1 

0 

2 

2 

0 

3 
0 
2 
0 
0 
0 
2 

4 
2 
0 
0 
2 
2 

5 
5 
2 
0 
0 
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0 

0 

0 

0 

***********************  S0URce  ******************** 

ELECTRONS 
ENERGY  15.0 

**********************  OPTIONS  ******************* 
PULSE-HEIGHT  1  1 
NBINE  34  USER 

14.99999  10.1  9.57  8.94  8.3  7.67  7.01  6.37  5.74  5.12  4.48  3.84 
3.2  2.56  1.91  1.28 

1.020  0.979  0.910  0.839  0.765  0.694  0.622  0.549  0.480 
0.408  0.336  0.263  0.193  0.125  0.051  0.049  0.00001  0.0 
ELECTRON-FLUX  1  1 
NBINE  19  USER 

14.99999  1.020  0.979  0.910  0.839  0.765  0.694  0.622  0.549  0.480 
0.408  0.336  0.263  0.193  0.125  0.051  0.049  0.00001  0.0 

DOME- SOURCE  0.  0.  0.0  1.25 

**********************  OPTIONS  ******************* 

CUTOFFS  0.1  0.01 
HISTORIES  10000000 
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Appendix  3 

CEASE  Telescope  -  MCNPX  Flat  Disk  Proton  Source  Runs 


MCNPX  Run  1.  Disc  source  located  at  z=0.457  cm  along  telescope  axis; 
disk  source  radius  =  0.025  cm; 

9  MeV  proton  source 
10,000  histories 

Partial  listing  of  output  file:  input  file  image,  pulse  height,  flux,  and  energy  deposition  for 
cells  67  (DFT)  and  88  (DBT)  _ 

1“  Test  for  protons  CEASE  Flight  Sensor  9  MeV  rad=0. 025cm 

2-  C 

3-  C  Proton 

4-  C  transport  -  disk  source,  9  MeV  proton  40  deg  slant  at  z=.4572 

GEOMETRY  FILE  DETAILS  SAME  AS  SHOWN  IN  [FIRST  YEAR  REPORT] APPENDIX  3  OMITTED  HERE  FOR 
BREVITY 

462-  C 

^63-  c  Transport  protons,  neutrons,  muons,  photons,  pions,  kaons, 

464-  c  deuterons,  alphas 

465- 

466-  mode  h 

467-  cut :h  l.e+8  0.05 

468-  c  mode  h  n  j  p  / 

469-  c  Source  definition,  proton  source  (par=9) ,  located  on  surface 

470-  c  #47  which  is  the  plane  at  z=0.4572,  centered  at  0, , 0.  ,0.4572, 

471-  c  with  radius  0.025  cm,  energy  9  MeV,40deg  slant  incidence  [along  z) 

472-  SDEF  sur=47  pos=0.  0.  .4572  rad=Dl  ERG=9 .  WGT=1.0  par =9 

473-  dir=0. 7660444 

474-  SI1  .025 

475-  VOL  45 J  0.144715  3J  4.501818  27J  0.1311943  21 J  0.1333961 

476-  2J  0.1129316  0.4679637  14J  19.61581  61.74832  5J 

477-  139.2479  139.2432  103.1126  103.1126  103.1126 

478-  103.1126  42CF 

479-  c  No.  of  histories 

480-  NFS  10000 

481-  C 

482-  C  Materials 

483-  C 

484-  C*  nracc 

485-  Ml  29000  -.3  28000  -.7 

486-  C  Aluminum 

487-  M2  13027  -1.0 

488-  C  Tungsten 

489-  M3  74000  -.95  29000  -.015  28000  -.035 

490-  C  Gold 

491-  M4  79197  -1.0 

492-  C  Stainless  Steel 

493-  M5  26000  -.71  29000  -.17  25055  -.065  28000  -.005  14000  -.05 

494-  C  Conductive  Silicone  Elastomer 

495-  M6  28000  -.377  47000  -.373  14000  -.0947 

496-  6000  -.0810  8016  -.0539  1001  -.0204 

497-  C  PMMA 

498-  M7  6000  -.59985  8016  -.31961  1001  -.080538 

499-  C  Silicon 

500-  MB  14000  -1.0 

501-  C  Copper 

502-  M9  29000  -1.0 

503-  C 

504-  C 

505-  C 

506-  c  maximum  proton  energy  [MeV)  required  for  cross  section  table 

507-  PHYS :h  10. 

508-  c  maximum  neutron  energyCMeV)  required  for  cross  section  table 

509-  PHYS:n  10. 

510-  c  maximum  muon  energy  (MeV)  required  for  cross  section  table 

511-  PHYS: j  10. 

512-  c  maximum  photon  energy  [MeV)  required  for  cross  section  table 

513-  PHYS:p  10. 

514-  c  maximum  pion  energy  [MeV)  required  for  cross  section  table 

515-  PHYS:/  10. 

516-  C 
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517- 

c 

Tallies 

518- 

c 

proton  pulse  height  tallies 

519- 

F8:h 

67 

520- 

E8  0 

l.e-5  .5  1.  1.5  2.  2.3  3.  3.5  4.  4. 

521- 

7.  7.5  8.  8.5  8.99999  9. 

522- 

F38 :h  88 

523- 

E38 

0  l.e-5  .5  1.  1.5  2.  2.3  3.  3.5  4.  4 

524- 

7.  7.5  8.  8.5  8.99999  9. 

525- 

c 

proton  energy  deposition  tallies 

526- 

*F18 

:h  67  88 

527- 

C 

proton  flux  tallies 

528- 

F44: 

h  67  88 

529- 

C 

c  proton  energy  flux  tallies 

530- 

*F104 :h  67  88 

531- 

C 

cell  importances  for  protons 

532- 

imp:h  1  165R  0  6R 

533- 

c 

cell  importances  for  photons 

534- 

c 

imp:p  1  165R  0  6R 

535- 

c 

cell  importances  for  neutrons 

536- 

c 

imp:n  1  165R  0  6R 

537- 

c 

cell  importances  for  muons 

538- 

c 

imp: |  1  165R  0  6R 

539- 

c 

cell  importances  for  pions 

540- 

c 

imp:/  1  165R  0  6R 

541- 

ltally 

8 

nps  =  10000 

tally 

type  8  pulse  height  distribution. 

particle (s):  proton 

cell  67 

energy 
0 . OOOOE+OO 
1 . 0000E-05 
5 . 0000E-01 
1. OOOOE+OO 
1 . 5000E+00 
2. OOOOE+OO 
2 . 3000E+00 
3. OOOOE+OO 
3 . 5000E+00 
4. OOOOE+OO 
4 . 5000E+00 
5. OOOOE+OO 
5 . 5000E+00 
6.  OOOOE+OO 
6 . 5000E+00 
7. OOOOE+OO 
7 . 5000E+00 
8. OOOOE+OO 
8 . 5000E+00 
9. OOOOE+OO 
9. OOOOE+OO 
total 


0 . OOOOOE+OO  0.0000 
0 . OOOOOE+OO  0.0000 
1 . 00000E-03  0.3161 
3 . 20000E-03  0.1765 
2.90000E-03  0.1854 
5 . 13000E-02  0.0430 
2 . 71000E-02  0.0599 
3 .  41000E-02  0.0532 
1 . 62000E-02  0.0779 
1.22000E-02  0.0900 
1 . 03000E-02  0.0980 
3 . 90000E-03  0.1598 
1 . 10000E-03  0.3013 
1 . OOOOOE-Q4  0.9999 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
1 . 63400E-01  0.0226 


tally  type  8*  energy  deposition 
particle(s):  proton 


cell  67 


4 . 15596E-01  0.0243 


units 


units 


cell  88 

2 . 61874E-02  0.1198 

tally  type  8  pulse  height  distribution.  units 

particle(s):  proton 

cell  88 


energy 

0. OOOOE+OO 

0. 00000E+00 

0.0000 

1.0000E-05 

0. OOOOOE+OO 

0.0000 

5 . 0000E-01 

2 . OOOOOE-04 

0.7070 

1. OOOOE+OO 

8.00000E-04 

0.3534 

1 . 5000E+00 

7. 00000E-04 

0.3778 

2. OOOOE+OO 

3. 00000E-04 

0.5773 

2 . 3000E+00 

5. 00000E-04 

0.4471 

3. OOOOE+OO 

1 . 40000E-03 

0.2671 

3 . 5000E+00 

1 . 20000E-03 

0.2885 

4. 0000E+00 

1 . 00000E-03 

0.3161 

4 . 5000E+00 

8. 00000E-04 

0.3534 

5. OOOOE+OO 

6. OOOOOE-Q4 

0.4081 

5 . 5000E+00 

5. OOOOOE-04 

0.4471 

6. OOOOE+OO 

4 . 00000E-04 

0.4999 

6 . 5000E+00 

1 . OOOOOE-04 

0.9999 

7. 0000E+00 

0. OOOOOE+OO 

0.0000 

7 . 5000E+00 

0. OOOOOE+OO 

0.0000 

number 


mev 


number 
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8 . OOOOE+OO 
8.5000E+00 
9, OOOOE+OO 
9. OOOOE+OO 
total 


0 . OOOOOE+OO  0.0000 
O.OOGOOE+OG  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
8 . 5000GE-03  0.1080 


tally  type  4  track  length  estimate  of  particle  flux, 
particle (s) :  proton 


units  l/cm**2 


volumes 


cell;  67  88 

7 . 61920E-03  4 . 23G17E-02 


cell 

67 

3 . 57828E-01  0.0241 

cell 

88 

2.20099E-03  0.1326 

tally  type  4*  track  length  estimate  of  energy  flux, 

par  tide  (s):  proton 

units  mev/em**2 

volumes 

cell;  67  88 

7 . 61920E-03  4 , 23017E-G2 

cell  67 


2 . 02689E+00  0.0268 

cell  88 

5 . 68473E-03  0.1553Test  for  protons  CEASE  Flight  Sensor  9  MeV 

rad=0 .025cm 
C 


MCNPX  Run  2.  Disc  source  located  at  z=0.457  cm  along  telescope  axis; 
disk  source  radius  =  0.050  cm; 

9  MeV  proton  source 
40,000  histories 

Partial  listing  of  output  file:  input  file  image,  pulse  height,  flux,  and  energy  deposition  for 
cells  67  (DFT)  and  88  (DBT) 


1- 

2- 

3- 

4- 

5- 
462“ 

463- 

464- 

465- 

466- 

467- 

468- 

469- 

470- 

471- 

472- 

473- 

474- 

475- 

476- 

477- 

478- 

479- 

480- 

481- 


Test  for  protons  CEASE  Flight  Sensor  9  MeV  rad=0 . 05cm 
C 

C  Proton 

C  transport  -  disk  source,  9  MeV  proton  40  deg  slant  at  z=.4572 

C 

C 

c  Transport  protons,  neutrons,  muons,  photons,  pious,  kaons, 

c  deuterons ,  alphas 

mode  h 

cut:h  l.e+8  0.05 

c  mode  hn  |  p  / 

c  Source  definition,  proton  source  <par=9),  located  on  surface 
c  #47  which  is  the  plane  at  z=0.4572,  centered  at  0. ,0. ,0.4572, 
c  with  radius  0.05  cm,  energy  9  MeV,40deg  slant  incidence  (along  z) 
SDEF  sur=47  pos=0.  0.  .4572  rad=Dl  ERG=9 .  WGT=1.0  par ^9 

dir-0. 7660444 

SI1  .05 

VOL  45 J  0.144715  3J  4.501818  27J  0.1311943  21J  0.1333961 
2J  0.1129316  0.4679637  14J  19.61581  61.74832  5J 

139.2479  139.2432  103.1126  103.1126  103.1126 

103.1126  42J 

c  No.  of  histories 
NPS  40000 
C 
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482- 

C  Materials 

483- 

C 

484- 

C 

Brass 

485- 

Ml 

29000  -.3  28000  -.7 

486- 

C 

Aluminum 

487- 

M2 

13027  -1.0 

488- 

C 

Tungsten 

489- 

M3 

74000  -.95  29000  -.015  28000  -.035 

490- 

C 

Gold 

491- 

M4 

79197  -1.0 

492- 

C 

Stainless  Steel 

493- 

M5 

26000  -.71  29000  -.17  25055  -.065  28000  -.005  14000  -.05 

494- 

C 

Conductive  Silicone  Elastomer 

495- 

M6 

28000  -.377  47000  -.373  14000  -.0947 

496- 

6000  -.0810  8016  -.0539  1001  -.0204 

497- 

C 

PMMA  - 

498- 

M7 

6000  -.59985  8016  -.31961  1001  -.080538 

499- 

C 

Silicon 

500- 

M8 

14000  -1.0 

501- 

C 

Copper 

502- 

M9 

29000  -1.0 

503- 

C 

504- 

C 

505- 

C 

506- 

C 

maximum  proton  energy (MeV)  required  for  cross  section  table 

507- 

PHYS 

:h  10. 

508- 

c 

maximum  neutron  energy  (MeV)  required  for  cross  section  table 

509- 

PHYS 

:n  10. 

510- 

C 

maximum  muon  energy  (MeV)  required  for  cross  section  table 

511- 

PHYS 

:|  10. 

512- 

C 

maximum  photon  energy  (MeV)  required  for  cross  section  table 

513- 

PHYS 

:p  10. 

514- 

C 

maximum  pion  energy  (MeV)  required  for  cross  section  table 

515- 

PHYS:/  10. 

516- 

c 

517- 

c 

Tallies 

518- 

c 

proton  pulse  height  tallies 

519- 

F8:h 

67 

520- 

E8  0 

l.e-5  .5  1.  1.5  2.  2.3  3.  3.5  4.  4.5  5.  5.5  6.  6.5 

521- 

7.  7.5  8.  8.5  8.99999  9. 

522- 

F38: 

h  88 

523- 

E38 

0  l.e-5  .5  1.  1.5  2.  2.3  3.  3.5  4.  4.5  5.  5.5  6.  6.5 

524- 

7.  7.5  8.  8.5  8.99999  9. 

525- 

C 

proton  energy  deposition  tallies 

526- 

*Fl8:h  67  88 

527- 

C 

proton  flux  tallies 

528- 

F44 : 

h  67  88 

529- 

C 

c  proton  energy  flux  tallies 

530- 

*F104:h  67  88 

531- 

c 

cell  importances  for  protons 

532- 

imp:h  1  165R  0  6R 

533- 

c 

cell  importances  for  photons 

534- 

c 

imp :p  1  165R  0  6R 

535- 

c 

cell  importances  for  neutrons 

536- 

c 

lmp:n  1  165R  0  6R 

537- 

c 

cell  importances  for  muons 

538- 

c 

imp: |  1  165R  0  6R 

539- 

c 

cell  importances  for  pions 

540- 

c 

imp:/  1  165R  0  6R 

541- 

1 tally 

8 

nps  =  40000 

tally  type  8  pulse  height  distribution.  units 

particle{s):  proton 

cell 

67 

energy 

0. 

OOOOE+OO 

0 . 00000E+00  0.0000 

1. 

0000E-05 

0 . 00000E+00  0.0000 

5. 

0000E-01 

5 . 75000E-04  0.2085 

1. 

0000E+00 

1 . 45000E-03  0.1312 

1. 

5000E+00 

1 . 90000E-03  0.1146 

2. 

OOOOE+OO 

1 . 49500E-02  0.0406 

2. 

3000E+00 

1.01750E-02  0.0493 

3. 

OOOOE+OO 

1 . 61250E-02  0.0391 

3. 

5000E+00 

7 . 27500E-03  0.0584 

4. 

OOOOE+OO 

5 . 50000E-03  0.0672 

4. 

5000E+00 

4 . 72500E-03  0.0726 

5. 

OOOOE+OO 

1 . 80000E-03  0.1177 

5. 

5000E+00 

2 . 00000E-04  0.3535 

6. 

0000E+00 

0 .  OOOOOE+OO  0.0000 

number 
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6 . 5000E+00 
7 . OOOOE+OO 
7.5G00E+00 
8. OOOOE+OO 
8.50G0E+00 
9. OOOOE+OO 
9. OOOOE+OO 
total 


0 . OOOOOE+OO  0.0000 
0 . OQGOGE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
0. OOOOOE+OO  0.0000 
2.50GG0E-05  1.0000 
0. OOOOOE+OO  0.0000 
6.47000E-02  0.0190 


tally  type  8*  energy  deposition 
particle (s) :  proton 

cell  67 


cell  88 


1.69238E-01  0.0202 


1.30129E-02  0.0846 

tally  type  8  pulse  height  distribution, 
par  tide  |s):  proton 

cell  88 


energy 

0. OOOOE+OO 

0. OOOOOE+OO 

0.0000 

1.0000E-05 

0. OOOOOE+OO 

0.0000 

5 . 0GG0E-01 

2 . 5000GE-G4 

0.3162 

1. 0000E+00 

3 . 0G000E-04 

0.2886 

1.5000E+00 

3.00000E-04 

0.2886 

2. OOOOE+OO 

5 .00000E-04 

0.2236 

2 . 3000E+00 

2.50000E-04 

0.3162 

3. OOOOE+OO 

9.25000E-04 

0.1643 

3.5000E+00 

4 .750002-04 

0.2294 

4. OOOOE+OO 

3 . OOOOOE-04 

0.2886 

4.5000E+00 

3 . 2500GE-04 

0.2773 

5. 0000E+00 

3.0G000E-G4 

0.2886 

5 . 5000E+00 

2.5G0G0E-04 

0.3162 

6. 0000E+00 

2. 000Q0E-04 

0.3535 

6.5000E+00 

7.50000E-05 

0.5773 

7. OOOOE+OO 

0. OOOOOE+OO 

0.0000 

7 . 5000E+00 

0. OOOOOE+OO 

0.0000 

8 . 00G0E+00 

0. OOOOOE+OO 

0.0000 

8 . 5000E+00 

0. OOOOOE+OO 

0.0000 

9. 0000E+00 

0. 00000E+00 

0.0000 

9. 0000E+00 

0. OOOOOE+OO 

0.0000 

total 

4.4500GE-03 

0.0748 

tally  type  4  track  length  estimate  of  particle  flux. 

particle {s};  proton 

volumes 

cell:  67  88 

7.61920E-03  4.23017E-02 

cell  67 

1.35987E-01  0.0203 

cell  88 


1.08395E-03  0.0956 

tally  type  4*  track  length  estimate  of  energy  flux. 

particle (s) :  proton 

volumes 


cell  67 


cell:  67  88 

7 . 61920E-03  4 . 23G17E-Q2 


cell  88 


7.09309E-01  0.0227 
2 . 81644E-03  0.1149 


units 


units 


units 
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Appendix  5 

MCNPX  Source  Subroutine  for  Grazing  Angle  Proton  Scattering  Study 
c_deck  so  source 
so  1 

subroutine  source 
so  2 

c  user  supplied  source  subroutine 
♦include  "cm.h" 
c 

c  This  is  the  source  routine  for  oblique  proton  beams 

c  incident  on  Al  or  Ir  slab 

c  The  name  of  this  deck  is  srcoblq.F 

c 

data  pin/3.14159265/ 
data  komin/0/ 
if (komin.eq. 0) then 
komin=l 

open { 47 ,  f  orm= ' unformatted ' , s tatus= ' scratch ' ) 
write ( jtty, 20) 

20  format (lx, 'Enter  source  energy  (MeV)  and  incident 
obliquity  (deg.) 

$') 

c  example,  2  degree  obliquity  means  enter  92.0  deg  (i.e. 
92  deg  off  normal) 

read (i tty, *) erg, obliq 
write (47) erg, obliq 
write (jtty, 21) erg, obliq 
write ( iuo ,21) erg , obliq 

21  format (lx, ' Source  energy  (MeV) =', el2 . 5 , 2x, ' Incident 
obliquity  (deg 

$.)  = ' , el2 .5) 
rewind  47 
end  if 
wgt=l . 0 
tme=0 . 0 

read ( 47 ) erg, obliq 
rewind  47 
c 

c  The  source  energy 

c  erg=  0.1 

c  ipt=9  denotes  proton  source, 

c 

ipt=9 

c  incident  grazing  obliquity  92  degrees  off  normal  (2 
degrees  off  surface) 

ths=obliq/180 .  *pin 
sths=sin(ths) 
cths=cos (ths) 

c  incident  azimuth  =90  degrees 

phs=0 . 5*pin 
sphs=l. 0 
cphs=0 . 0 
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c  direction  cosines  of  incident  beam 

uuu= s  ths  *  cphs 
vw= s  ths  *  sphs 
www=cths 

aa=sqrt  (uuu**2+vw**2+www**2) 

uuu=uuu/aa 

wv=vw/aa 

www=www/ aa 

c  point  of  source  incidence 

xxx=0 . 0 
yyy= -0.99 
z z z =yyy * c ths / s t hs 
jsu=0 
icl=2 

do  50  ispr=l,3 

50  spare (ispr) =0 . 0 

return 
end 
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Appendix  6 

PHILOOK,  Program  for  Extraction  and  Conversion  to  Histogram  Form  of 
Emergent  Proton  Azimuthal  Distributions  from  MCNPX  PTRAC  File  Output 
Produced  in  the  Grazing  Angle  Proton  Scattering  Study 


implicit  real*8  (a-h,o-z) 

dimension  phibin(60) , count (60) ,a{200) ,b(200) 
character* 120  line 
data  pye/3. 14159265/ 
data  count/60*0./ 

open(2,file='lookit' , status = 'unknown' ) 
open(7,file='binl' , status= ' unknown ' ) 
write (6, 1) 

1  format (lx, 'Enter  number  of  phi  bins  and  number  of 
histories ' ) 

read ( 5 , * ) nbins , his ts 
dphi=360 . /nbins 
write (6,17) nbins 
17  format (lx, i5) 

phibin ( 1 ) =dphi 
do  10  n=2 , nbins 

10  phibin (n) =phibin(n-l) +dphi 

fac=180 .  /pye 

open(l, file='ptrac' , status= ' old' ) 
do  20  i=l, 10 
20  read(l,2) line 

2  format (a) 

100  read(l, 2 , end=200) line 

readd,  2)  line 

readd,  *)x,y,  z,u,  v,  w,  e,wt,  t 

s th=sqrt ( 1 . -w*w) 

if (sth.eq. 0 . ) go  to  100 

cphi=u/sth 

sphi=v/sth 

if (cphi.ge. 0.  .and.  sphi . ge. 0 . ) kq=l 
if (cphi.ge.0.  .and.  sphi . It . 0 . ) kq=4 
if (cphi . It . 0 .  .and.  sphi . ge . 0 . ) kq=2 
if (cphi.lt. 0.  .and.  sphi.lt.0. )kq=3 
if (abs (sphi-1 . 0) .le.l.e-5)kq=5 
if (kq.eq.5)phi=0 . 5*pye 
if (kq.eq.l)phi=asin(abs (sphi) ) 
if (kq.eq.2)phi=pye-asin(abs (sphi) ) 
if  (kq.eq.3)phi=pye+asin(abs  (sphi) ) 
if (kq.eq.4)phi=2 . *pye-asin(abs (sphi) ) 
phid=fac*phi 
do  30  n=l , nbins 


if (phid.gt .phibin(n) )go  to  30 
ibin=n 

if (ibin.eq.l)write(7,77)x,y, z,u, v,w, sphi , cphi 
77  format (lx, 8el6 . 9) 

count ( ibin) =count ( ibin) +1 . 
go  to  100 
30  continue 

200  do  210  n=l,nbins 

count (n) =count (n) /hists 
210  write(2/4)n,phibin(n) , count (n) 

4  format (lx, 'bin  no . = ' , i5, lx, 'upper 
bound (deg. )  = ' , f7 .2 , lx, ' count= ' ,  e 

$12.5) 

i=l 

a(l)=0. 

b(l) =count (1) 

nbl=nbins-l 

do  215  n=l , nbl 

i=i+l 

a(i)  =phibin(n) 
b(i) =count  (n) 
i=i+l 

a(i)=phibin(n) 

215  b(i) =count (n+1) 
ii=i+l 
a (ii) =360 . 
b(ii) =b(l) 
do  220  i=l,ii 
220  write(2,5)a(i) ,b(i) 

5  format (lx, f7 . 1,  2x,el2.5) 
close  (1) 

close  (2) 

stop 

end 
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